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Abstract 

We calculate the 3-loop perturbative expansion of the average pla- 
quette in lattice QCD with Nf massive Wilson fermions and gauge 
group SU(N). The corrections to asymptotic scaling in the corre- 
sponding energy scheme are also evaluated. We have also improved 
the accuracy of the already known pure gluonic results at 2 and 3 
loops. 
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1 Introduction and calculation 



The lattice numerical calculation of the matrix element of any operator A 
with mass dimension d yields a d A, where a is the lattice spacing, up to 
scaling violating contributions of order 0(a d+1 ). The calculation of a as a 
function of the coupling constant can be improved if an effective scheme is 
used 0. A nice example of such improvement has been shown to occur in 
pure Yang-Mills theories || and is known to work rather well in other scalar 
theories (see for example [J| for the 2D spin models). 

In this paper we calculate the 3-loop internal energy in full QCD with 
Wilson fermions on the lattice. The 3-loop beta function for this theory 
has been recently calculated |], 0] and we can compute the corresponding 
corrections to asymptotic scaling in both the standard and effective schemes. 

We will use the Wilson action S with Nf flavours and gauge group SU (N) . 
This action is 



S 
Sw 



Sw + Sf, 



with 

E f (x) 



l--ReTr(D) 



E - 

flavours 



( m + ~a) ^ x ^ : 



(2) 



(3 is the coupling on the lattice (3 = 2N/g 2 and r the Wilson parameter. 
□ stands for the plaquette and U^(x) = exp (iagA^x)) for the link at the 
site x pointing towards x + jl. We will parametrize the results in terms of N, 
Nf and the couple (k, r) where k is the usual hopping parameter 
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8r + 2am 



(3) 
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The average of Ef can be straightforwardly computed to all orders by 
rescaling the fermionic action Sf by a factor e under the fermionic path 
integral Z* 

Z f (e) = J Vi>(x)Vi)(x) exp (-eS f ) = e iVNN fZ f (e = 1) (4) 
and by using 

m = 4 i^) ri - -*™,. (5) 

On the other hand, the average of E w in presence of the action Eq. @ is 
calculated in perturbation theory 

(Ew) = c 1 g 2 + c 2 g i + c 3 g 6 + --- (6) 

The n-loop coefficient can be written as c n = c 9 n + c{ where c 9 n is the pure 
Yang-Mills contribution, known since ref. || Q up to 3-loops, and c{ is the 
fermionic contribution. To calculate c{ we will first compute the free energy 
— (\nZ)/V up to 3 loops, Z being the full partition function 

Z = J VU^{x)Vi,(x)V4)(x)exv(-S). (7) 

The average of Ew is then extracted as follows 

The Feynman diagrams necessary to calculate the fermion contribution 
to the free energy up to 3 loops are shown in Fig. 1. We have used the 
Feynman gauge. The involved algebra of the lattice perturbation theory was 
carried out by making use of the computer code developed by us 

After grouping the diagrams in several infrared-finite sets, we calculated 
the resulting finite integrals on finite lattice-sizes L and then the results 
were extrapolated to infinite size. The extrapolating function was of the 
type i@ 

ao+ 2^ Oij-j^-. (9) 

i<2j, j=l,2,- u 
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We used a broad spectrum of such functional forms and analyzed the quality 
of each extrapolation to assign a weight to each one, to finally produce a 
reliable estimate of the systematic error. The criterion to judge the quality of 
the extrapolation was based on the accuracy of the the fitted functional form 
to reproduce known results at finite lattice-sizes. The different functional 
forms were obtained by truncating the series in Eq.flj^) at different values of 
j and assuming vanishing coefficients for some i and j. 

Recall || |9| that the pure gluonic contributions are (presented here with 
improved accuracy) 



r 9 
c l 



r 9 

c 2 



r 9 
c 3 



N 2 -l 
8N 



(iV 2 - l) (o.0051069297 - — L_) , (10) 

/ 2 \ /0.0023152583(50) 0.002265487(17) 
V ~ ) { JP N + 

0.000794223(19) AM. 



The 2-loop coefficient can be written in closed form as 

Pi 2 Qi Q 2 i 



(iV 2 




12 24 288 128 A^ 2 



where Pi, Q\ and Q2 are finite integrals defined and evaluated in | |13| |. 

In tables I and II we show c 2 and c{, for several pairs {K,r). They are 
parametrized in terms of four constants h%, h 30 , h^i and h 32 as follows 



c{ = 



J 

c 2 



(iV 2 - l) h 2 (12) 



4 = (N 2 - l) (h m N f + h 31 ^ 



N 2 




In table III we show the result for c 2 = cf + c 2 and C3 = cf + c{ for r = 1, 
N f = 3 and A^ = 2. In table IV the result for r = 1, N f = 3 and N = 3 is 
shown. 
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2 Corrections to asymptotic scaling 

The beta function in QCD with Nf Wilson fermions can be written as 

t{g) = , = ~bo9 3 - hg 5 - b 2 g 7 - • • • (13) 

The non-universal 3-loop coefficient b 2 has been recently calculated in ref. |J, 
(see also 0). In terms of the bare coupling g, the lattice spacing a approaches 
the continuum limit as 

aA L = exp (Vr i/26 ° (l + 19 2 + • ■ ■) , (14) 



where g is the 3-loop correction to asymptotic scaling 

_ bj - b 2 b 

q = —tiT- (15) 

Other couplings can be defined. A popular effective coupling in terms of 
the plaquette energy is [|l|, 

9e w = ^ (Ew)mc, (16) 

where (-)mc indicates Monte Carlo average. In terms of qe w the approach of 
the lattice spacing to the continuum limit is written as in Eq.([TJ|) with q^ w 
instead of q. This 3-loop correction to asymptotic scaling is 

b c 3 - bi c 2 - b cl/ct 
2 Cl bl 

Moreover the lattice Lambda parameters are related by the equation 



QE w = q • (I 7 ) 



A ^ = exp (2^) Ai ' (18) 

We give the coefficients for SU(2) in table V and for SU(3) in ta- 
ble VI. In both cases we show the result for r = 1 and for several choices of 
k. These results must be compared with the value of q in the standard bare 
coupling scheme. This is [[J 

SU(2) -> q(N f = 1) = 0.12617(4) q(N f = 3) = 0.2547(1), 
S77(3) -»■ g(JV> = 1) = 0.23956(4) q(N f = 3) = 0.3681(2), (19) 
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(notice that in the standard bare coupling scheme, q does not depend on the 
fermion mass 0). The 3- loop correction to asymptotic scaling for Nf = 3, 
N = 3 is ~ 37% in the standard scheme and ~ 18 — 21% in the (E w ) 
scheme. Apparently this improvement is not as dramatic as it was for the 
pure gluonic case || (for SU(3) q ~ 19% in the standard scheme and q ~ 1% 
in the energy scheme). The improvement therefore seems to be more efficient 
in the quenched case. This fact can also be seen from tables V and VI 
where the approach to the quenched case (either Nf — > or k — > 0) is 
accompained by the lowering of qE w - However, the only relevant test for 
the improvement in asymptotic scaling from the standard to some energy 
scheme is the practical use of the effective scheme. For example, the 2D spin 
models [^|, [J show a definite improvement in spite of the behaviour of the 
3-loop coefficient q which for the 0(3) models passes from q = —0.09138 in 
the standard scheme to 0.1694 in the energy scheme. 

3 Discussion 

We have calculated the free energy up to 3 loops in QCD with Wilson 
fermions. We have given the result of the internal energy average (Ew) 
as a function of the number of fermions Nf, their masses and the Wilson 
parameter r. These expansions have been used to study the 3-loop correc- 
tions to asymptotic scaling in this theory. We have found that at 3 loops 
the corresponding energy scheme |l], |2|] provides a moderate improvement 
with respect to the standard scheme. However, only the practical use of this 
scheme in particular problems will reveal how useful it is. 

We can numerically compute the set of coefficients and q also for other 
choices of (k, r), Nf and N. 
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Figure caption 



Figure 1. Feynman diagrams contributing to the fermionic sector of InZ/V at 2 
and 3 loops. Curly, arrowed and dashed lines indicate gluons, fermions 
and ghosts respectively. 



Table captions 

Table I. 2-loop fermion contribution to the internal energy (E w ). See text for 
notation. 

Table II. 3-loop fermion contribution to the internal energy (E w ). See text for 
notation. 

Table III. Coefficients C2 and C3 for SU(2) and Nf — 3 at r = 1. 

Table IV. Coefficients C2 and C3 for SU(3) and Nf = 3 at r = 1. 

Table V. 3-loop correction to asymptotic scaling in SU (2) with r = 1 in the 
energy scheme (Ew)- 

Table VI. 3-loop correction to asymptotic scaling in SU (3) with r = 1 in the 
energy scheme (Ew)- 
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Table I 



K 


r 


h 2 x 10 3 


0.1675 


0.1 


-0.269079(27) 


0.161 


0.1 


-0.236105(26) 


0.156 


0.1 


-0.212480(25) 


0.1675 


0.5 


-0.502547(23) 


0.161 


0.5 


-0.434010(22) 


0.156 


0.5 


-0.385777(21) 


0.1675 


1 


-2.408379(14) 


0.164 


1 


-2.273850(12) 


0.16 


1 


-2.117310(36) 


0.1575 


1 


-2.017993(24) 


0.156 


1 


-1.957882(14) 


0.1675 


2 


-6.44013(90) 


0.161 


2 


-6.29318(41) 


0.156 


2 


-6.16472(95) 



Table II 



K 


r 


h 30 x 10 3 


h 31 x 10 3 


h 32 x 10 3 


0.1675 


0.1 


-0.0194886(93) 


0.030311(54) 


0.000585(20) 


0.161 


0.1 


-0.017144(10) 


0.026733(47) 


0.000451(18) 


0.156 


0.1 j 


-0.015466(11) 


0.024140(41) 


0.000366(17) 


0.1675 


0.5 


-0.032691(24) 


0.050837(49) 


0.002041(15) 


0.161 


0.5 


-0.028269(21) 


0.044149(52) 


0.001521(16) 


0.156 


0.5 


-0.025205(15) 


0.039464(51) 


0.001202(14) 


0.1675 


1 


-0.361568(64) 


0.4271834(82) 


0.049019(10) 


0.164 


1 


-0.330133(54) 


0.390773(14) 


0.043856(15) 


0.16 


1 


-0.293659(46) 


0.348654(20) 


0.038177(14) 


0.1575 


1 


-0.270560(43) 


0.322111(22) 


0.034772(11) 


0.156 


1 


-0.256620(43) 


0.306132(22) 


0.0327854(66) 


0.1675 


2 


-1.4249(32) 


1.6754(19) 


0.34483(95) 


0.161 


2 


-1.3838(44) 


1.6257(34) 


0.3297(14) 


0.156 


2 


-1.3541(29) 


1.5868(22) 


0.3166(14) 
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Table III 



K 


c 2 x 10 3 


c 3 x 10 3 


0.1675 


-1.376291(62) 


0.60413(60) 


0.164 


-0.770912(55) 


0.73543(54) 


0.16 


-0.06648(16) 


0.89226(47) 


0.1575 


0.38045(11) 


0.99446(44) 


0.156 


0.650945(61) 


1.05715(42) 




Table IV 






K 


c 2 x 10 3 


c 3 x 10 3 




0.1675 


14.64396(11) 


7.3440(16) 




0.164 


15.720192(98) 


7.8775(14) 




0.16 


16.97251(29) 


8.5042(12) 




0.1575 


17.76705(19) 


8.9061(12) 




0.156 


18.24794(11) 


9.1504(12) 
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Table V 



K 


q Ew (N f = l) 


q Ew (N f = 3) 


0.1675 


0.05318(5) 


0.2052(2) 


0.164 


0.05069(4) 


0.1954(1) 


0.16 


0.04782(5) 


0.1842(2) 


0.1575 


0.04599(5) 


0.1770(1) 


0.156 


0.04489(5) 


0.1727(1) 


Table VI 


K 


q Ew (N f = l) 


QE w (N f = 3) 


0.1675 


0.06645(4) 


0.2091(2) 


0.164 


0.06355(4) 


0.1989(2) 


0.16 


0.06021(4) 


0.1871(1) 


0.1575 


0.05809(4) 


0.1796(1) 


0.156 


0.05681(4) 


0.1752(2) 
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Fig. la 



